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BOUNDING  GLOBAL  MINIMA  WITH 
INTERVAL  ARITHMETIC 

by 

* 

Louis  J.  Mancini 
Garth  P.  McCormick 


It  is  shown  how  techniques  of  interval  arithmetic  can  be  used  to 
give  bounds  on  the  global  minima  of  unconstrained  optimization  problems. 
The  techniques  are  illustrated  using  the  design  of  a hypothetical  chemical 
plant. 
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1.  Introduction 

Many  problems  in  all  branches  of  engineering,  business  and  economics 
can  be  mathematically  formulated  as  a nonlinear  programming  problem  (NLP) . 

A NLP  is  an  optimization  problem  in  which  some  objective  is  minimized  (or 
maximized)  subject  to  certain  constraints  which  describe  the  system  being 
modeled.  The  objective  may  be  to  minimize  costs,  as  the  example  in  this 
section  shows,  or  to  maximize  some  performance  level.  The  adjective 
"nonlinear"  refers  to  the  fact  that  the  objective  and  constraint  functions 
may  be  nonlinear  in  contrast  to  the  functions  used  in  linear  programming. 

In  this  respect,  linear  programming  can  be  viewed  as  a special  case  of  non- 
linear programming;  although  the  solution  procedures  used  in  each  are 
considerably  different. 

If  no  constraints  are  present,  then  the  NLP  is  termed  "unconstrained" 
and  can  be  symbolically  written  as 


min  f(x) 
x 


(1) 
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where  f , the  objective  function,  is  a scalar  function  of  n variables 
(i.e.,  the  vector  x ) . For  example,  consider  the  design  of  the  hypo- 
thetical chemical  plant  shown  in  Figure  1.  The  design  variables  are 

x^  ■ reactor  temperature 

x2  * reactor  pressure 

x^  * weight  fraction  of  a catalyst  used 

x^  * weight  fraction  of  the  product 
leaving  the  reactor. 

The  product  x^  is  an  explicit  function  of  the  temperature  x^  given  as 
0.25 

x^  “ 0.1  x^*  . Summing  the  component  costs  in  Table  1 and  substituting 

for  x^  yields  the  cost,  or  objective,  function 

f(x)  - 0.0318  xj*1  x®*6  + 11430  x'1  x"1  + 228  x3  + 1.5  x2  - 25  xj*25  - 495  xj’2  . 


Table  1 

COMPONENT  COSTS  FOR  CHEMICAL  PLANT 


Item 

Cost 

($/100  lb.  material 
processed) 

1.  Reactor 

0.0318  xj’1  x®*6 

2.  Separator 

11430  x'1  X31 

3.  Catalyst 

228  x3 

4.  Compressor 

1.5  x2 

5.  Recycle  compressor 

250  (1  - x4) 

6.  Byproduct  sales 

<)|  AA  A 0.8 
-3123.2  x. 

4 

- 2 - 


Figure 


T-345 


NLPs  are  usually  solved  with  an  algorithm  which  starts  with  a given 
point  (i.e.,  the  best  guess)  and  performs  some  iterative  procedure  until 
the  solution,  or  an  approximate  solution,  is  found.  A difficult  problem 
in  any  such  algorithm  is  termination,  that  is,  when  should  the  iterative 
procedure  be  stopped?  Computer  programs  which  implement  algorithms  must 
contain  some  "convergence  criterion"  which  is  testable  and  which  if  satis- 
fied causes  termination  of  the  algorithm.  When  this  happens  the  algorithm 
is  said  to  have  "converged." 

One  class  of  convergence  criteria  are  those  which  use  the  iterates 

0 k k+1  , . . , , 

x ,...,x  ,x  ....  to  determine  whether  or  not  convergence  has  been 

achieved.  The  user  is  required  to  supply  a value  C (presumed  small) 

before  the  algorithm  is  initiated.  The  algorithm  is  terminated  when 

I k+1  k I r 

x - x < 


or 


|f(xk+1)-f(xk) | < C 


depending  upon  the  criterion  selected.  The  obvious  flaw  in  the  use  of 
these  criteria  is  that  closeness  between  the  successive  iterates  does  not 
imply  the  iterate  is  close  to  x*  , the  solution  to  (1).  The  main  advantage 
of  these  criteria  is  that  they  prevent  bad  algorithms  (ones  which  fail  to 
change  the  iterates  very  much  each  iteration  or  fail  to  decrease  the  objective 
function  significantly  each  time)  from  using  a lot  of  computer  time. 

From  the  user’s  point  of  view,  it  is  desired  to  terminate  only 
when  the  iterate  or  its  function  value  is  acceptably  close  to  x*  (or  f(x*)) 
That  is,  from  considerations  of  the  physical  meaning  of  the  problem  it  is 
reasonable  to  ask  the  user  to  supply  a value  C and  for  the  algorithm  to 
terminate  when  either 

i k 


**l 


< G 


or 


|f(xk)-f(x*)|  < 
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In  some  rare  cases  the  value  f(x*)  is  known  even  though  x*  is 
not.  Thus  the  fourth  criterion  would  be  implementable  in  this  circumstance. 

In  general  though,  this  is  not  the  case  and  these  natural  convergence  cri- 
teria cannot  be  used.  It  is  a major  contribution  of  this  paper  that  using 
the  theoretical  material  developed  in  [7]  with  the  tools  of  interval  arith- 
metic [5],  that  these  criteria  can  be  implemented  for  unconstrained  NLPs. 
Before  showing  this  it  is  instructive  to  consider  another  convergence  cri- 
terion which  has  been  suggested  in  an  attempt  to  bypass  the  lack  of  know- 
ledge of  x*  . 

A necessary  condition  for  a vector  x*  to  be  the  optimal  solution 
to  an  unconstrained  NLP  is  that 

f'(x*)  = £ (2) 

where  f'  is  the  gradient  vector  and  £ is  the  zero  vector.  Suppose  we 
are  solving  (1)  with  some  iterative  procedure  and  the  algorithm  has  gen- 

Jq 

erated  a vector  x . Based  on  (2),  we  could  terminate  if  the  norm  of  the 

i k i 

gradient  vector  is  near  zero,  that  is,  if  |f'(x  ) | <_ C . In  terms  of 

k k k k 

the  chemical  plant  example  suppose  x = (x^,  x^,  x^)  = (698.33,  35.738, 

1.1929)  and  C=  0.05.  Then  |f'(xk)|  = 3.25  and  we  would  not  terminate. 
However,  now  suppose  we  measure  costs  in  dollars/lb  rather  than  dollars/100  lb 
(see  Table  1),  and  thus  define  c(x)  = f(x)/100  as  the  new  cost  function. 

Then  for  the  same  vector  xk  , |c'(xk)|  = |f'(xk)|/100  = 0.0325  , and  we 
would  terminate.  Thus  this  criterion  fails  since  the  norm  of  the  gradient 
vector  is  dependent  upon  the  scaling  of  the  problem. 

There  are  two  further  objections  to  the  use  of  this  last  convergence 
criterion.  Whereas  the  user  may  have  a good  reason  to  accept  a point  which 
is  within  C in  norm  from  the  solution  vector  or  whose  objective  function 
is  within  C of  the  optimal  objective  function  value,  there  is  usually  no 
physical  interpretation  of  a point  whose  gradient  norm  is  less  than  C . 
Second,  as  a practical  matter  no  one  ever  tries  this  convergence  criterion 
more  than  once.  Except  for  simple  problems,  because  of  numerical  considera- 
tions, there  Heems  to  be  a value  above  zero  below  which  it  is  impossible  to 
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reduce  the  gradient  norm  even  with  the  most  powerful  minimization  techniques. 
The  computer  just  keeps  running  and  running. 

This  paper  presents  a solution  to  the  termination  problem  for 
unconstrained  NLPs.  In  [7]  the  authors  have  shown  that,  under  certain 
conditions,  an  analytic  expression  exists  for  the  difference  f(x)  - f(x*)  ; 
where  x*  is  the  global  solution  over  some  compact  region.  This  difference 
represents  how  far,  in  terms  of  the  objective  function,  the  vector  x is 
from  the  optimal  vector  x*  . However,  the  necessary  conditions  and  the 
analytic  expression  for  the  difference  usually  cannot  be  computed  exactly. 

It  will  be  shown  how  interval  analysis  can  be  used,  under  certain  conditions, 
to  compute  an  interval  bound  on  the  difference  f(x)  - f(x*)  . In  the  pro- 
cess an  interval  bound  on  the  solution  vector  x*  is  also  generated.  These 
bounds  can  be  used  as  termination  criteria. 

2.  Interval  Analysis 

Interval  arithmetic  [8]  generalizes  ordinary  arithmetic  to  closed 
intervals  of  the  real  line.  Give  two  intervals  u = [a,b]  and  w = [c,d]  , 
interval  arithmetic  is  defined  by 

u + w = [a  + c,  b + d] 

u - w = [a  - d,  b - c] 

u • w = [min(ac,  ad,  be,  bd) , 
max(ac,  ad,  be,  bd)] 

u/w  = [a,  b]  • [1/d,  1/c] 

where  u/w  is  defined  only  if  0 i [c,d]  . The  degenerate  interval, 

[u,u]  , is  not  distinguished  from  the  ordinary  number  u . On  the  computer, 
rounded  interval  arithmetic  [6,15]  is  used  in  place  of  the  exact  version 
above  to  bound  round  off  error.  Rounded  intervals  might  be  slightly  larger, 
but  they  are  guaranteed  to  contain  the  exact  result. 

Interval  matrices  are  rectangular  arrays  with  intervals  as  components. 
A square  interval  matrix  A is  nonsingular  if  and  only  if  all  ordinary 
matrices  AcA  are  nonsingular,  that  is, 

- 6 - 
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0 i {det  A|AeA}  . 

Usually  it  is  not  possible  to  compute  the  right-hand  side  of  the  above 
exactly,  however,  an  Interval  determinant,  det  A,  satisfying 

{det  A|AeA}  C Tet  A (3) 

can  be  computed.  If  A is  nonsingular,  then  an  interval  inverse,  (A) 
satisfies 

{A-1|AeA}  C (A)'1  . 

Consider  the  interval  linear  system 

A y = b (4) 

where  A is  nonsingular.  The  theoretical  solution  set,  SI  , is 

ft  = {y|Ay  = b,  AeA,  beb) 

and  an  interval  solution  y must  satisfy  b C a y ; and  hence  contains  fi  . 

The  y with  smallest  width  [8,  p.  7]  is  desired  to  minimize  the  difference 
between  A y and  b . Hansen  {3,4],  Hansen  and  Smith  [5,6],  and  Oettli  [10] 
have  developed  methods  for  solving  (4). 

If  f:Rn-*R  is  continuous,  then  an  interval  extension  of  f , on  an 

interval  vector  x C Rn  , is  an  interval-valued  function  f which  satisfies 

f(x)  * f(x)  for  all  xex 

(5) 

f(u)  C f (w)  for  all  u C w C x 

Interval  extensions,  like  interval  determinants  and  interval  inverses,  are 
not  unique;  and 

{f(u)|ueu}  C f(u)  (6) 

for  all  u Q,  x . The  "best"  interval  extension,  or  united  extension  [8,  p.  18], 
is  that  function  f in  which  equality  holds  in  (6).  Finding  good  computable 
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interval  extensions,  limiting  the  difference  between  the  left-  and  right- 
hand  sides  of  (6),  is  one  of  the  key  problems  in  interval  analysis. 

Suppose  h:Rn-»-Rn  is  continuously  differentiable,  and  let  x*  denote 
a zero  of  h . Assume  continuous  [8,  p.  18]  interval  extensions,  , of 

the  partial  derivatives 


hIj(x)  = 9x~  ’ i,j  = 

are  available  on  x . For  u,wC_  x define 

h^(u,w)  = h^j(u1,...,uj_1,wj,...,wn)  , i,j  = 1, . . . ,n  (7) 

and  let  h' (u,w)  be  the  interval  matrix  with  components  h'.  . (u ,w)  . If 

h'(x,x)  is  nonsingular  in  the  interval  sense,  then  h'(u,w)  is  non- 
singular for  all  u,wC  x ; and  an  interval  Newton  operator,  N^  , can  be 
defined  as 

N^w.w)  = w - h'(w,w)  1 • h(w)  for  all  w e w C x . (8) 

Nickel  has  proven  that  the  above  interval  Newton  operator  has  some  interesting 
properties. 

Theorem  1 [9] : Suppose  h'(x,x)  is  nonsingular.  Then 

(a)  Any  zero  x*  of  h in  x is  unique. 


(b)  If  x*  e x , then  x*  e N^(x,x)  for  all  x e x . 

(c)  If  xO  N^x.x)  = 0 for  any  x e x , then  x*  t x . 

(d)  If  N^(x,x)  C x for  any  x e x , then  x*  e N^x.x)  . 


If  the  hypotheses  in  Part  (d)  hold,  then  the  existence  of  a zero  x*  of  h 
in  x has  been  verified,  x*  is  contained  in  the  interval  N^x.x)  , and 
x*  is  the  unique  zero  of  h in  x . 
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13.  Bounding  Global  Minima 

The  authors  have  shown  [7]  that,  under  certain  conditions,  an 
analytic  expression  exists  for  the  difference  f(x)  - f(x*)  ; where 

* 

f:Rn+R  and  x*  minimizes  f over  a convex  compact  set,  S , in  Rn  . 
However,  the  necessary  conditions  and  the  expression  for  the  difference 
usually  cannot  be  computed  exactly.  This  section  presents  a computational 
procedure,  based  on  interval  analysis,  which  can  produce  an  interval  bound 
on  the  difference  f(x)  - f(x*)  . In  the  process  an  interval  bound  on 
x*  is  also  generated. 

Assume  f is  twice  continuously  differentiable  on  S , and  that 
the  Hessian  matrix,  f’’(x)  , is  positive  definite  on  S . Define  the 
gradient  set 

U = {u|u  = f'(x)  for  some  x e S}  . 


Ml 


T-345 


(iii)  e S‘  , or 
(iv)  N(x°)  C S1  ; 

then  tf ' (x^)  e U for  all  t e [0,1]  , and 

f(x°)  - f(x*)  = f’(xV  • /J  tf”[g{tf’(x0)}]_1  dt  • f’(x°)  . (9) 

At  first  glance  the  above  theorem  seems  to  be  of  only  theoretical  interest. 
Hypotheses  (i),  (ii)  and  the  expression  (9)  for  the  difference  f(x°)  - f(x*) 
usually  cannot  be  computed  exactly.  However,  it  will  be  shown  that  interval 
analysis  can  be  used  to  calculate  an  interval  version  of  Theorem  2. 

Suppose  the  convex  compact  set  S is  some  interval  x in  R , 
and  assume  continuous  interval  extensions,  f|!  , of  the  second  partial 

derivatives  of  f are  available  on  x . Let  f''(x,x)  be  the  interval 
Hessian  matrix  with  components  f!J(x,x)  given  by  (7).  If  an  interval 

determinant  (3)  is  computed,  and 

0 t dit  f' r(x,x)  , (10) 

then  f”(x)  is  nonsingular  for  all  x e x . Furthermore,  since  the 
eigenvalues  of  a matrix  are  continuous  functions  of  the  matrix  components 
[11];  if  f"(x)  is  positive  definite  for  any  x c x and  (10)  holds, 
then  f"(x)  is  positive  definite  on  x , and  hypothesis  (i)  of  Theorem  2 
is  satisfied. 

An  interval  Newton  operator  can  be  used  to  verify  hypothesis  (ii)  . 

The  definition  of  an  interval  extension  (5)  implies  that 

/*  f"[x°  + (x  - x°)s]ds  e f''(x,x) 

for  all  x°  , x in  x . If  (10)  holds,  then  f''(x,x)  is  invertible; 
and  the  definition  of  an  interval  inverse  implies  that 
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(/J  f”[x°  + (x-x°)s]ds)-1  e f'Cx.x)"1  (11) 

for  all  x^  , x In  x . An  interval  Newton  operator,  N2  , which  uses  the 
symmetric  f''(x,x)  can  be  defined  as 

N2(x°,x)  = x^  - f"(x,x)  1 * f'(x°)  . (12) 

Therefore,  if  x^  e x and 

N2(x^,x)  C_  x (13) 

then  (11)  and  (12)  imply  that 

N(x°)  C_N2(x°,x)  C x , 
and  hypothesis  (ii)  is  satisfied. 

If  (10)  and  (13)  hold,  then  from  Theorem  2,  there  exists  a stationary 

point  x*  of  f in  N2(x^,x)  , and  x*  is  the  unique  such  point  in  x . 

Furthermore,  if  f''(x)  is  positive  definite  for  some  x c x , then  x* 
minimizes  f over  x . Assume  now  that  either  the  point  x®  and/or 

the  interval  N2(x^,x)  lie  in  the  interior  of  x . Then  from  Theorem  2, 

tf'(x^)eU  for  all  tc[0,l]  , where  U = {u|u  *=  f'(x)  for  some  x e x } . 
Therefore,  the  definitions  of  an  interval  extension  and  interval  inverse 
imply  that 

f'fgttfMx0)}]  e f"(x,x) 

and 

/J  tf’MgUfHx0)}]"^  eif"(x,x)-1  . 

Using  the  above  with  (9)  yields  an  interval  bound  on  the  difference  in 
function  values 

f(x°)  - f(x*)  e \ f'(x°)T  • f”(x,x)_1  • f'(x°)  . (14) 


- 11  - 
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The  previous  discussion  has  proven  the  following  interval  version 
of  Theorem  2. 

Corollary;  Let  x be  some  interval  in  Rn  , let  x*  denote  the 
interior  of  x , and  assume  continuous  interval  extensions,  f^  , of  the 

second  partial  derivatives  of  f are  available  on  x . Suppose  x^*  e x 
and 


i! 


f 


i 


( i ' ) 0 i det  f"(x,x)  and  f''(x)  is  positive 

definite  for  some  x e x , 

( ii f ) the  interval  Newton  operator  (12)  is  contained 
in  x , that  is,  N2(x^,x)  C x . 

Then  there  exists  a x*  e N2(x^,x)  such  that 

f ’ (x*)  - g , f(x*)  = inf^  f(x)  , 

xex 

and  x*  is  the  unique  stationary  point  of  f in  x . Furthermore,  if 
(iii')  x^  e x1  , or 
(iv')  N2(x^,x)C  x1  ; then 

f(x°)  - f(x*)  e 1/2  f'(x°)T  • f”(x,x)  f'(x°)  . 

The  right-hand  side  of  the  above  is  an  interval  bound  on  the  difference 

f(x^)  - f(x  ) , and  the  interval  Newton  operator  (12)  is  an  interval  bound 
on  x*  . 

Note  that  the  first  part  of  the  corollary  pertaining  to  the  existence 
and  uniqueness  of  x*  is  similar  to  Parts  (a)  and  (d)  of  Theorem  1 by 
Nickel  (let  h(x)  ■ f'(x)  in  Theorem  1).  Actually,  Part  (d)  of  Theorem  1 
is  stronger  since  it  uses  a tighter  interval  Newton  operator.  This  follows 
from  the  definition  of  an  interval  extension  since  (5)  and  (7)  imply  that 

7*  • (x°,x)  c 7’  • (7.7) 
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for  all  x°  e x , and  since  interval  arithmetic  is  inclusion  monotonic 
[8,  p.  10] 

N^(x°,x)  C N2(x°,x) 


where 


N1(x°,x)  - x°  - f”(x°,x)  1 • f'(x°)  . 


Therefore,  the  interval  Newton  operator  N^(x^,x)  , which  uses  the  non- 


symmetric  f"(x°,x)  , might  be  contained  in  x in  cases  where  N2(xw,x) 
defined  by  (12)  is  not. 

Consider  the  computations  involved  in  implementing  the  corollary. 

Given  x in  Rn  , x°  t x , and  continuous  interval  extensions,  fJJ  , 

on  x ; the  symmetric  f''(x,x)  is  formed,  and  det  f''(x,x)  is  calculated. 
An  interval  determinant,  det  A , can  be  calculated  by  diagonalizing  A , 
using  partial  pivoting  in  interval  arithmetic,  to  obtain  an  upper  triangular 


rr  , 0 -v 


A ; then 
u 


det  A - x (A  )..  . 
i=l  U 


It  is  easy  to  show  that  the  above  satisfies  (3).  If  (10)  holds,  then 
f"(x)  is  nonsingular  on  x , and  the  interval  Newton  operator  N2(x°,x) 

given  by  (12)  is  well  defined.  Now  N2(x°,x)  can  be  calculated  by  first 

solving  the  interval  linear  system 

f ' ' (x,x)  • y * -f ' (x°) 

for  an  interval  solution  y , and  augmenting  x^  to  obtain 

N2(x^,x)  » x°  + y . (16) 


The  procedure  used  here  to  solve  for  y is  Hansen  and  Smith's  Method  4 [5]. 
In  [4]  Hansen  presents  a refinement  procedure  which  may  improve  a given 
interval  solution. 
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If  (13)  holds,  then  ^(x^.x)  contains  a stationary  point  x*  of  f , 

and  x*  is  the  unique  such  point  in  x . Furthermore,  if  f''(x)  is  posi- 
tive definite  for  any  x e x , then  f''(x)  is  positive  definite  on  x ; 
and  x*  minimizes  f over  x . The  interval  bound  (14)  on  the  difference 

f(x°)  - f(x*)  holds  if  x^*  e x1  and/or  ^(x^.x)  C x1  , and  is  easily 

obtained  from  the  calculation  of  N2(x^,x)  . This  follows  since  (12),  (14), 

and  (16)  imply  that 

x°  - N2(x°,x)  = f(x,x)-1  • f'(x°)  = - y 

and  therefore 

f(x°)  - f(x*)  e - f'(x°)T  • y . 

It  cannot  be  overemphasized  how  much  the  usefulness  of  the  above 
results  depends  on  the  availability  of  good  interval  extensions  and  good 
numerical  methods  for  calculating  an  interval  determinant  and  solving  an 
interval  linear  system. 


Unconstrained  Signomial  Programs 


Many  optimization  problems,  especially  in  engineering  design  [13], 
can  be  formulated  as  an  unconstrained  nonlinear  program  where  the  objective 
function  f , a signomial  [12],  can  be  written  as 


f(x) 


(18) 


Here  x is  a positive  vector  in  Rn  , and  the  coefficients  cf  and 

exponents  afci  are  arbitrary  real  numbers.  If  all  the  coefficients  c^ 

are  positive,  then  f is  a posynomial  [2];  and  has  a unique  stationary  point 

which  is  the  global  minlmizer.  Defining  the  exponent  vectors 

a^  ■ (flti , . . . »atn^  an<*  substituting  z • log  x , f can  be  rewritten  as 

m 

f(z)  - E s (z) 
t-1 
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I 

I 


where  the  single-term  functions,  s^  , are  defined  as 

T 

8t(«)  “ c exp(at  • z)  , t=l,...,m  . 

For  an  interval  [u,w],  define  exp([u,w])  » [exp(u) ,exp(w) ] ; and 
for  positive  [u,w] , define  log([u,w])  * [log(u) ,log(w) ] . 

Given  a positive  vector  x , the  corollary  requires  an  interval 
extension  of  the  Hessian 

f"(z)  = " s (z)  a • a l 

t*l  c t t 

on  z = log  x . Since  each  s^  is  monotonic,  an  interval  extension 

st(z)  , is  easily  obtained  by  computing  in  interval  arithmetic.  Then 

the  definition  of  an  interval  extension,  the  inclusion  monotonicity  of 
interval  arithmetic,  and  Theorems  4.1  and  4.2  by  Moore  [8,  p.  19]  imply 
that 

m _ _ 

f”(z,z)  = l s (z)  a,.  • a^  (19) 

t-1  c C C 

is  a continuous  interval  extension  of  f’’(z)  evaluated  on  z . The  above 
procedure  is  computationally  attractive,  however,  different  methods  exist 
which  might  yield  sharper  results. 

As  an  example  consider  the  minimum-cost  design  of  the  hypothetical 
chemical  plant  introduced  in  Section  1.  The  procedure  will  be  illustrated 
on  two  intervals  z * log  x containing  the  (approximate)  stationary  point 


(694.9\  /6.S437  \ 

35.56  ) - ( 3.5712  ) (20) 

1.187/  \0. 17143/ 

obtained  using  Dembo's  algorithm  [1]  for  signomial  programs.  This  will  permit 
a comparison  between  the  interval  bounds  on  the  difference  f(z^)  - f(z*) 
with  the  actual  values.  Note  that  since  f''(z*)  is  positive  definite, 
if  f''(z,z)  is  nonsingular  in  the  interval  sense,  then  f''(z)  is  positive 
definite  on  z . 
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Table  2 presents  results  where  the  interval  on  the  original  design 
variables  x equals  a + 2.5  percent  symmetric  deviation  around  x*  given 
in  (20)  , that  is, 

([677.48,  712. 24]  \ 

[34.671,  36.450]  ) . 

[1.1573,  1.2167]/ 

Substituting  z = log  x and  computing  the  interval  Hessian  (19)  yields 

([338.20,  376.31] 

[229.05,  249.39] 

0 

Calculating  the  interval  determinant  defined  by  (15)  yields 

dTt  7" (1,1)  = [1.6802(107),  4.6121(107)]  . 

Therefore,  f''(z)  is  positive  definite  on  z , and  the  interval  Newton 
operator 

N2(z°,7)  - z°  - f"(7J)-1  • f'(z°)  (21) 

is  well  defined.  As  discussed  in  the  previous  section,  N2(z°,z)  is 

calculated  by  solving  the  interval  linear  system 

f"(z,z)  • y « - f’(z°)  (22) 

for  an  interval  solution  y and  augmenting  z®  to  obtain  N^z^.z)  = z^*  + y . 

The  results  for  three  different  points  z^  = log  x^  (where  x^  = x*  + 0.005  x*  , 
x*  + 0.01  x*  , x*  + 0.02  x*)  are  given  in  Table  2,  and  in  each  case 

N2(z<7,z)Cz  . Therefore,  in  each  case  the  existence  of  a stationary  point 

3 

z*  in  N2(z  ,z)  has  been  verified,  z*  is  the  unique  stationary  point  in 

z and  z*  minimizes  f over  z . Furthermore,  since  in  each  case  z^ 
lies  in  the  interior  of  z , the  interval  bound 

f(z°)  - f(z*)  c - | f’(z°)T  • y (23) 

holds;  and  since  z*  is  known  in  this  example,  the  bound  can  be  compared 
with  the  actual  value.  Note  that  as  the  distance  between  z®  and  z* 


[229.05,  249.39]  0 

[434.68,  475.56]  [257.74,  284.86] 

[257.74,  284.86]  [521.61,  562.26] 
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becomes  larger,  the  interval  solution  y to  (22)  becomes  wider;  ana  hence 
the  interval  Newton  operator  and  the  interval  bound  become  wider. 

Table  3 presents  results  with  the  same  points  x^  and  where  the 
interval  x equals  a + 5 percent  deviation  around  x*  , that  is, 


([660.11,  729.61]  \ 

[33.782,  37.338]  I . 

[1.1276,  1.2464]  / 


Substituting  z = log  x and  computing  (19)  yields 


f " (z,z) 


([319.64,  395.86] 

[219.15,  259.81] 

0 


[219.15,  259.81] 
[415.82,  497.78] 
[245.61,  300.06] 


[245.61,  300.06]  ] . 

[502.72,  584.23]  / 


Again,  f''(z)  is  positive  definite  on  z since  (15)  yields 


det  f"(z,z)  = [0.41092(107) , 6.3177(10?)]  . 

However,  the  requirement:  N2(z^,z)  C z only  holds  in  the  first  case  where 
z^  = log  (x*  + 0.005  x*)  . In  cases  2 and  3,  the  combination  of  the  width 
of  z and  the  distance  between  z^  and  z*  cause  the  method  to  fail. 

All  basic  interval  arithmetic  operations  were  performed  using  Zoltan's 
single-precision  rounded  interval  arithmetic  package  [15].  Other  numerical 
procedures  yielding  tighter  interval  extensions,  or  better  methods  for 
solving  an  interval  linear  system  should  produce  even  better  results. 
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